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INTRODUCTION 


The  radar  and  computer  system  can  be  considered  as  the  eyes 
and  brain  of  a  ballistic  missile  defense  (BMD)  system.  Given  no  limi¬ 
tation  on  resources,  a  ballistic  missile  system's  operational  philoso¬ 
phy  may  be  to  utilize  predetermined  battle  plans  for  the  allocation  of 
radar  and  computer  resources.  However,  it  is  rarely  economically 
feasible  to  provide  excessive  resources,  and,  since  the  defensive  end 
offensive  situation  changes  dynamically,  it  becomes  imperative  that 
the  allocation  of  such  tesources  adapt  to  these  changes. 


This  paper  treats  the  radar  and  computer  Bystem  as  an  adap¬ 
tive  control  system  in  which  the  allocation  of  the  computer  preprocess¬ 
ing  rate,  radar  pulse  rate,  and  computer  postprocessing  rate  are 
varied  dynamically.  To  determine  the  optimal  real-time  allocation  of 
the  radar  and  computer  resources,  the  problem  was  first  formulated 
mathematically  and  then  solved  in  two  phases. 


PROBLEM  FORMULATION 


To  develop  quantitative  as  well  as  qualitative  results,  a 
typical  BMD  system  problem  is  considered.  The  principal  physical 
elements  of  the  problem  are  a  radar,  a  computer,  interceptors,  and 
threatening  reentry  vehicles  with  decoys.  The  radar  and  the  inter- 
captor  battery  are  under  the  control  of  a  computer.  The  problem  is 
then  to  allocate  the  two  major  resources,  i.e.,  the  radar  and  the  com¬ 
puter,  to  the  task  of  searching  for  potential  future  threats,  tracking 
identified  objects,  launching  Interceptors,  guiding  Interceptors, 
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tracking  interceptors,  and  processing  all  related  data.  A  typical 
objective  of  this  allocation  is  to  maximize  the  probability  that  the 
radar  survives  the  attack.  Any  similar  criteria,  subject  to  the 
constraints  on  the  available  resources,  such  as  minimizing  the  proba¬ 
bility  of  a  reentry  vehicle  penetrating  the  defet  3e,  may  also  be  con¬ 
sidered  as  the  objective. 

The  measurement  of  the  achievement  of  the  objeitive  of  the 
problem  can  be  considered  a  function  of  the  following  variables: 

•  The  number  of  radar  pulses  allocated  to  searching 

•  The  number  of  objects  presently  detected  and  the  possi¬ 
bility  of  detecting  future  objects 

•  The  number  of  radar  tracking  pulses  allocated  to  the  1th 
object 

•  The  impact  points  of  the  threatening  objects 

•  The  possibility  of  an  object  being  a  decoy 

•  The  reliability  cf  the  interceptor 

•  The  possible  intercept  points  of  the  engagement 

•  The  number  of  interceptors  allocated  to  each  target* 

•  The  number  of  radar  guidance  pulses  allocated  to  the  jth 
interceptor  which  has  been  committed  to  the  ith  object 

•  Natural  and  man-made  environmental  effects. 

Based  on  the  sample  problem,  l.e.,  to  maximize  the  probability  that 
the  radar  survives  thii  attack,  the  control  variables  of  the  problem 
are  as  follows: 

•  The  search  rate  tor  future  objects,  F 

•  The  future  track  rates  of  the  n  objects  presently 
detected,  Ki,  ...,  K„ 

•  The  number  of  track  and  guidance  pulses  for  each  inter¬ 
ceptor  assigned  to  the  n  objects,  L},  ...,  Ln 

The  major  constraints  of  the  sample  problem  are  the  available 
radar  power,  the  computer  processing  speed,  and  the  scheduling  of 
Kj,  ...,  Kn;  Lj,  •••*  ^n* 

The  sample  problem  can  then  be  stated  as  follows: 

Maxluize:  (0)  Probability  the  radar  survives  ■  f  (F,  K^,  Lj) 

n  n 

Subject  to:  (1)  auF  +  au  <2  Kj  +  al3  (2  L.)  s  b»(t) 

i-1  *  i-1  1 


*Each  target  is  assigned  a  single  interceptor. 
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n  n 

(2)  a2iF  +  a22  (2  K.)  +  823  (2  L.)  <  b2(t) 

i-1  1  i-1 

(3)  Scheduling  of  Ki,  ....  Kn;  Lj,  Lj, 

subject  to  time  frame  limitations, 

where  bi(t)  -  Maximum  amount  of  radar  energy  available  per  second 

b2(t)  “  Maximum  number  of  basic  computer  instructions  available 
per  second 

an  *  Radar  energy  required  for  a  single  search  pulse 

an  *  Radar  energy  required  for  a  single  target  radar  track 

pulse 

au  ■  Radar  energy  required  for  a  single  Interceptor  radar 
pulse 

a2i  “  Number  of  basic  computer  Instructions  required  to 
process  each  search  pulse 

a22  “  Number  of  basic  computer  instructions  required  to  pro¬ 

cess  each  target  track  pulse 

a23  “  Number  of  basic  computer  instructions  required  to  pro¬ 

cess  each  interceptor  track  pulse 

The  probability  of  radar  survival  (Equation  0)  must  be  expressed  quan¬ 
titatively  as  a  function  of  the  search,  discrimination,  and  tracking 
algorithms,  and  other  associated  functions.  Thus,  given  a  particular 
HMD  system's  characteristics,  one  can  decompose  the  probability  of 
survival  into  various  subfunctions  which  can  be  determined  analyti¬ 
cally. 

The  equations  of  the  first  two  constraints  are  formed 
straightforwardly;  however,  the  third  constraint,  i.e.,  scheduling, 
bns  a  nonlinear,  time-varying  nature  and  cannot  easily  be  expressed  in 
an  analytical  form.  Therefore,  it  would  be  difficult  to  solve  the 
problem  directly  by  considering  all  three  constraints  simultaneously, 
but  it  is  possible  to  separate  the  problem  into  two  phases.  In  Thase 
I,  the  allocation  problem  is  solved  without  considering  the  third  con¬ 
straint,  i.e.,  assuming  that  the  transmitted  and  receive  radar  and 
computer  times  are  contiguous.  In  Phase  II,  the  required  computer  and 
radar  times  that  were- determined  in  Phase  I  are  scheduled  such  that 
the  overall  schedule  time  is  minimized. 


PHASE  I 

Phase  I  is  used  to  establish  the  optimal  blend  of  resources 
for  the  problem  formulation  which,  as  stated  previously,  is: 
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Maximize:  (0)  Probability  the  radar  survives  -  f  (F,  K^,  L^) 

n  n 

Subject  to:  (1)  qiF  +  a32  (2  K.)  +  au  (2  L  )  s  bi(t) 

i-1  1  i-1  1 

n  n 

(2)  a21F  +  a22  (2  K  )  +  a23  (2  L  )  «  b2(t) 
i-1  i-1  1 

F,  Kif  L±  *  0 

and  is  solved  using  the  nonlinear,  constrained  search  method  denoted 
as  the  Complex  method  [1]. 

The  following  example  problem  with  n  threatening  objects  was 
solved  by  the  Complex  method  programmed  in  FORTRAN  and  executed  on  a 
Control  Data  Corporation  7600  computer  system.  Although  the  problem 
is  somewhat  smaller  than  would  be  encountered  in  practice,  it  does 
exhibit  the  typical  form  of  real-life  problems. 

The  performance  criterion  for  measuring  the  achievement  of 
radar  survivability  aan  be  expressed  as: 

n 

Z  =  PD(F)  {R  Pjd^)} 

i-1 

where 

Pd(F)  =  Probability  of  detecting  all  future  threats,  as  a  func¬ 
tion  of  search  rate  F 

PiptKj)  *  Probability  of  survival  against  the  ith  object  attacking 
the  radar,  as  a  function  of  object  track  rate 

Pl(Li)  =  Probability  of  survival  against  the  ith  object  attacking 
the  radar,  as  a  function  of  track  rate  of  the  inter¬ 
ceptor  assigned  to  ith  object 

n  *  Number  of  threatening  objects 

The  assumed  typical  performance  functions  are: 

Pd(F)  -  1.0  -  0.8  exp(-F/1.5) 
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RjO^)  =  0.95  -  0.7  exp(-Ki/2.?) 

PI(Li)  «  1.0  -  0.6  exp(-Li/3.5) 

These  functions  are  both  monotonic  (more  resource  produces  more 
performance)  and  concave  (more  resource  means  that  Incremental  per¬ 
formance  gain  decreases)  and  therefore,  it  can  be  shown  that  the 
necessary  and  sufficiency  conditions  for  a  unique  global  optimum  are 
satisfied.  The  problem,  with  assumed  values  for  radar  energy  per 
pulse,  computer  instructions  per  pulse,  and  total  radar  and  computer 
resources,  becomes: 

Maximize:  Z  -  {l.0-0.8  exp(-F/1.5)(  (0.95-0.7  exp(-K/2.5)|n 

)l. 0-0.6  exp(-L/3.5)}n 

Subject  to:  0.1F  +  n(0.02K)  +  n(0.003L)  <1.0 
1.0F  +  n(0.2K)  +  n(0.3L)  <  20.0 
F,  K,  L  2  0 
n  ■  1,  2,  « « « ,  10 

The  solution  obtained  is  given  in  Table  1,  along  with  the  computer 
execution  time  for  each  solution. 

PHASE  II 


Phase  II  is  used  to  establish  the  optimal  schedule  of  tne  al¬ 
location  results  obtained  in  Phase  I  for  the  radar  pulse  rate  and  com¬ 
puter  processing  rate  so  as  to  minimize  total  system  processing  time. 
The  problem  of  optimally  scheduling  tasks  (n  pulse  group, i)  through  the 
computer  and  radar  is  treated  as  a  "two '  machine  problem  with  feedback 
That  is,  tasks  are  first  passed  through  the  computer  (machine  1)  for 
preprocessing,  to  the  radar  (machine  2)  for  transmission  and  reception, 
and  then  back  to  the  computer  for  postprocessing. 

This  scheduling  problem  is  a  special  case  of  the  flow-shop 
sequencing  problem  which  has  been  the  topic  of  considerable  research 
in  the  past  few  years  [2] [3] [A] [5]  6] [7].  However,  because  of  the  com¬ 
binatorial  nature  of  the  problem  (nt  combinations)  which  results  in 
large  storage  requirements  and/or  long  computational  times,  most  work 
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has  been  limited  to  small  problems  or  to  heuristic  methodology.  To 
Illustrate  the  Immense  size  and  nature  of  the  scheduling  problem,  con¬ 
sider  a  20-job  problem.  If  one  could  evaluate  the  solution  to  each 
permutation  (there  are  201)  In  a  picosecond,  It  would  take  more  tha.: 

28  days  to  try  all  possibilities. 

Another  solution  method  using  a  mixed  Integer  formulation 
modified  for  evaluation  with  the  branch-and-bound  procedure  his  been 
Investigated  [7].  This  method  produces  a  linear  programming  problem, 
with  from  3n  to  a  maximum  possible  3(2n-l)  constraints  and  with  3n+l 
variables,  that  must  be  evaluated  In  order  to  obtain  a  solution  at  each 
branch-and-bound  node.  It  is  clear  that  this  method  Is  not  suited  for 
large  problems  because  of  both  storage  requirements  and  computational 
time  requirements.  Other  methods  and  examples  of  the  scheduling 
problem  can  be  found  In  Reference  [8]. 

The  two-machine  problem  with  feedback  Is  a  special  case  of 
the  flow-shop  problem,  but  It  retains  enough  similarity  that  heuristic 
scheduling  appears  to  be  the  most  feasible  method  of  approach.  Two 
heuristic  algorithms  have  been  developed  especially  for  solving  this 
scheduling  problem  and  they  hold  for  the  following  assumptions: 

e  All  jobs  and  machines  are  simultaneously  available. 

e  All  jobs  are  of  equal  importance. 

e  No  jobs  can  be  processed  simultaneously  by  more  than  one 
machine. 

•  No  machine  can  simultaneously  process  more  than  one  job. 

•  Process  times  are  independent  and  deterministic. 

e  Transportation  time  Is  either  included  in  the  process 
time  or  Is  negligible. 

Now,  define  the  process  time  of  job  a  on  machine  m  to  be  ta  K 
for  the  first  round  of  job  a  and  t’a>m  for  the  second  roun<T.  Also,  * 
consider  a  partial  schedule  o_  (which  may  be  null)  and  a  job  a  which, 
when  augmented  to  o_,  yields  the  schedule  (perhaps  IncompleteT  oa.  Now, 
let  T(oa,m)  be  the  completion  time  of  partial  schedule  oa  at  machine 
m  (first  romd)  and  let  T'  (oa,m)  be  the  completion  time  of  the  partial 
schedule  (second  round)  on  machine  5.  (The  completion  time  of  the 
last  job  on  the  last  machine  is  often  called  the  MAKESPAN.) 

The  first  heuristic  algorithm  (Method  1)  developed  for  solv¬ 
ing  this  scheduling  problem  will  now  be  described. 

Step  1.  For  each  job  1,  calculate  f(i)  -  t^  ^/t^  ^ 


46 


SATTERFIELD 


Step  2.  Arrange  the  jobs  in  alternating  order  of  f(i;, 
beginning  with  the  smallest  f(i).  Let  the  com¬ 
plete  schedule  n  so  generated  be 

S  -  (a;,  B2,  ....  an). 

Step  3.  Calculate  the  completion  times  T(a. ,1)  and 
T'(ai,2). 

Step  A.  If  T(a  ,1)  +  t'  .  >  T(a  ,2),  indicating 
ii  f  n 

that  the  schedule  obtained  is  optimal ,  go  to 
Step  11;  otherwise,  continue  with  Step  5. 

n 

Step  5.  Calculate  A  ■  £  (t  ,  +  t'  . ) 

i.l  ai*!  “i*1 

and  B  ■  £  /t  ,\  +  min  It  .  +  t'  .1 

Step  6.  If  A  >  B,  go  to  Step  11;  otherwise,  continue 
with  Step  7. 

Step  7.  Find  the  smallest  t.  ,  and  the  smallest  t' 

ai.l  aJ  »1 

If  i  f  J,  set  job  a^  first  in  the  schedule  and  a. 
last,  and  then  go  to  Step  10;  otherwise,  continue 
with  Step  8. 

Step  8.  Find  the  smallest  t.  i  (k  +  i)  and  the  smallest 

k» 

t  af ,1  (f  j  j). 

Step  9.  If  t  -  +  t'  .  <  t  -  +  t*  .,  set  job  a, 

*i,l  af»l  ak»x  aj»!  1 

first  in  tl.j  schedule  and  a^  last;  otherwise,  set 
job  a^  first  in  the  schedule  and  aj  last,  and 
continue  with  Step  10. 

Step  10.  Arrange  the  remaining  jobs  in  alternating  order  of 
f(l),  beginning  the  second  scheduled  job  with  the 
smallest  f(l).  Let  the  final  schedule  n  so  gener¬ 
ated  be 

S  ■  (a j ,  a2,  • • • ,  aQ) . 


kj 


SATTERFIELD 


Step  11.  Calculate  the  MAKESPAN  of  the  schedule  so 
obtained. 


The  second  heuristic  algorithm  (Method  2)  developed  for  solv¬ 
ing  this  scheduling  problem  Is  now  described*. 


Step  1.  For  each  job  1,  calculate 


f(D 


(ti .i  -  y> , 

min  (t  ”,  t  .) 
i(i  i»- 


Step  2.  Arrange  tha  jobs  In  ascending  order  of  f(i).  Let 
the  complete  schedule  n  so  generated  be 

S  ■  (aj,  a2,  •••»  *Q) • 

Step  3.  Calculate  the  completion  times  T(a. ,1)  and 
T(«i»2) . 

Step  4.  If  T(a  ,1)  -ft'  >  T(a  ,2),  Indicating  that 
n  *1, i  n 

the  schedule  Is  optimal,  go  to  Step  6;  otherwise, 
continue  with  Step  5. 


Step  5.  Find  1  such  that  T(a^,2)  -  T(aQ,l).  Arrange  jobs 
1  +  1  through  n  In  ascending  order  of  t*  2*  Let 
the  final  schedula  n,  as  obtained,  be  * 


S  ■  (a},  a2,  •••»  • 

Continue  with  Step  6. 


Step  6.  Calculate  tha  MAKESPAN  of  the  schedule  so  obtained. 

Repeat  Steps  4  and  5  until  no  further  improvement 
is  possible  or  until  an  optimal  solution  Is 
obtained. 


^Development  of this  algorithm  was  based  on  suggestions  offered  by 
Dr.  J«  N.  D.  Gupta,. 
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ALGORITHM  COMPARISON 

A  branch-and-bound  algorithm  has  been  programmed  for  compar¬ 
ing  the  effectiveness  of  the  two  heuristic  scheduling  algorithms.  As 
stated  earlier,  the  branch-and-bound  method  requires  large  storage  and 
long  computational  time.  However,  it  was  felt  that  for  the  class  of 
problems  that  it  could  handle  some  inference  could  be  made  as  to  the 
efficiency  of  the  heuristic  methods. 

The  branch-and-bound  method  developed  to  solve  this  schedul¬ 
ing  problem  is  an  adaptation  of  one  by  Lomnicki  [2].  Each  node  in  the 
branch-and-bound  method  repzesents  a  sequence  of  from  1  to  n  jobs. 
Consider  node  P,  corresponding  to  the  partial  schedule  or  where  or 
contains  a  particular  subset  (of  size  r)  of  the  n  jobs.  Then  a  lower 
bound  on  the  MAKESPAN  of  all  schedules  that  begin  with  sequence  o  is: 

'’•i.i)  ] 

where  T'(or,l)  and  T(cr,2)  are  the  completion  times  for  the  last  job 
of  the  r  jobs  in  the  sequence  on  machine  1  (second  round)  and  machine 
2,  respectively,  and  or  is  the  set  of  n-r  jobs  that  have  not  been  ; 
assigned  a  position  in  the  sequence. 

The  LB  is  first  evaluated  for  the  n  classes  of  permutations, 
i.e.,  for  those  starting  with  1,  2,  ...,  n,  respectively.  From  the 
vertex  with  the  lowest  value,  the  LB  is  evaluated  for  the  n-1  sub¬ 
classes.  From  these  values  and  the  previously  calculated  values,  the 
lowest  vertex  value  is  found,  from  which  the  LB  is  again  evaluated. 
Proceeding  in  this  manner  and  after  a  finite  number  of  steps,  a 
MAKESPAN  LB*  is  found  such  that  all  the  remaining  vertices  of  the 
scheduling  tree  have  values  greater  than  LB*.  The  permutation  with 
this  label  LB*  gives  an  optimal  solution  to  the  scheduling  problem. 

Some  preliminary  results  will  now  be  discussed.  A  total  of 
69  problems  have  been  solved  by  all  three  methods.  The  problems  were 
all  constructed  with  machine  times  (tj^)  selected  randomly  from  a 
uniform  distribution  of  numbers  from  1  to  50  for  machine  1  and  from  1 
to  100  for  machine  2. 

First,  a  comparison  of  the  average  computational  times  for 
the  three  methods  is  shown  in  Figure  1.  As  expected,  the  branch-and- 
bound  method  required  more  time  than  either  of  the  heuristics  and  is 


LB(or)  -  MAX 


+i3r  ( 


T(or,2) 
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not  consistent  with  real-time  control  requirements,  increasing  rapidly 
as  the  number  of  jobs  n  increases.  The  highest  average  time  for  either 
heuristic  was  less  than  0  1  second  for  n  ■  200,  which  is  within  real¬ 
time  control  constraints. 

Next,  heuristic  Method  1  found  the  optimal  solution  in  66 
percent  of  the  cases,  whereas  heuristic  Method  2  found  the  optimal 
solution  in  45  percent  of  the  cases.  Heuristic  Method  1  found 
solutions  equal  to  or  lower  than  Mithod  2  in  85  percent  of  the  problems 
solved.  However,  in  no  case  did  the  solutions  by  either  heuristic 
method  exceed  the  optimal  solution  by  more  than  3  percent.  More  work 
is  needed  to  evaluate  the  effects  of  the  distributions  and  ranges  of 
machine  times  on  the  heuristic  methods'  efficiency  and  selection. 

CONCLUSION 

A  two-phase  method  for  solving  the  resource  allocation  and 
scheduling  problem  of  an  adaptive  control  system  such  as  that  of  a  HMD 
system  has  been  developed  and  is  under  evaluation.  The  allocation  and 
scheduling  have  been  placed  into  an  analytical  context  from  which  it 
is  possible  to  derive  a  dynamic  methodology  that  can  allocate  BMD 
resources  in  an  effective  manner. 

In  Phase  I,  e  nonlinear,  constrained,  search  method  wa::  shown 
to  solve  a  typical  example  of  a  BMD  resource  allocation  problem  involv¬ 
ing  search,  interceptor  tracking,  and  object  tracking  in  less  than  0.1 
second.  This  solution  time  is  well  within  the  constraints  for  a  real¬ 
time  control  BMD  system. 

In  Phase  II,  the  scheduling  problem  was  structured  as  a  two- 
machine  problem  with  feedback  so  as  to  minimize  total  system  process¬ 
ing  time.  Because  of  the  often  immense  size  of  the  problem,  heuristic 
scheduling  was  employed,  with  two  competing  algorithms  developed.  A 
preliminary  investigation  of  the  efficiency  of  the  algorithms  was  per¬ 
formed  by  using  a  branch-and-bound  algorithm  for  finding  the  solution 
to  restricted  cases.  The  investigation  indicates  that  the  heuristic 
method  solutions  and  computational  times  are  consistent  with  real- 
time  control  requirements. 

It  is  concluded  that  the  scheme  that  has  been  developed, 
which  gives  an  analytical  approach  for  structuring  and  relating 
resource  functions,  ie  adequate  for  real-time  performance  criteria. 
Besults  to-date  indicate  the  practicality  of  utilizing  this  scheme 
under  the  stringent  real-time  requirements  of  a  typical  BMD  system. 
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